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In this paper we study isotropic integrable systems based on the braid-monoid algebra. 
These systems constitute a large family of rational multistate vertex models and are 
realized in terms of the B n , C n and D n Lie algebra and by the superalgebra Osp(n\2m). 
We present a unified formulation of the quantum inverse scattering method for many of 
these lattice models. The appropriate fundamental commutation rules are found, allowing 
us to construct the eigenvectors and the eigenvalues of the transfer matrix associated to 
theB n , C n , D n , Osp(2n-l\2), Osp(2\2n-2), Osp(2n-2\2) andOsp(l\2n) models. The 
corresponding Bethe Ansatz equations can be formulated in terms of the root structure 
of the underlying algebra. 
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1 Introduction 



In this paper we look at the problem of diagonalization of the transfer matrix of a certain 
class of integrable two-dimensional lattice models. Their Boltzmann weights are intimately 
connected with a rational Baxterization of the braid-monoid algebra (see e.g. ref. [|lj for a 
review) . In particular, we are going to analyse multistate vertex models which are based on 
the symmetries B n , C n , D n , and Osp(n\2m). 

One possible method of finding the eigenvalues of a given transfer matrix is by using 
the so-called analytical Bethe Ansatz 0. This technique relies on the unitarity, crossing 
and analyticity properties of the transfer matrix and, in some cases, an extra amount of 
phenomenological input is also required . This method has been applied to some of the models 
which we are going to consider in this paper, more precisely for the systems B n , C n , D n 
[D, £| and Osp(l\2n) || ||. Unfortunately, the explicit construction of eigenvectors of the 
transfer matrix is beyond the scope of the analytical Bethe Ansatz. The construction of exact 
eigenvectors, besides being an interesting problem on its own, is certainly an important step 
in the program of solving integrable systems. Thus, another route has to be taken if one wants 
to benefit from the knowledge of the eigenvectors. This would be the case of computing lattice 
correlation functions in the framework developed by Izergin and Korepin || |9|. 

A more powerful mathematical method, based on first principles, is the quantum inverse 
scattering method 0. This technique, together with the Yang-Baxter relation, offers us a 
unified viewpoint for studying the properties of integrability of two-dimensional solvable models 
[0, [H], TTH . One important feature of this method is that it permits us to present an 
algebraic formulation of the Bethe states. This step, however, depends much on our ability 
to disentangle the Yang-Baxter algebra in terms of appropriate commutation relations. The 
simplest structure of commutation rules has been discovered in the context of the 6— vertex 



model Tttt and its multi-state generalizations [12, 13 



In this paper, we shall deal with the diagonalization problem of the transfer matrix of certain 
rational braid-monoid vertex models by means of the quantum inverse scattering method. We 
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shall see that this program can be carried out in a universal way for a quite general class of 
systems: the B n , C n , Djft, Osp{2n - 1|2), Osp{2\2n - 2), Osp{2n - 2|2) and Osp{l\2n) vertex 
models. We show that the fundamental commutation rules have a common form in terms of 
the corresponding Boltzmann weights. As a consequence, the derivation of the eigenvectors, 
the eigenvalues and the associated Bethe Ansatz equations also have a quite general character 
for these vertex models. We believe that the unified picture proposed in this paper is a new 
result in the literature as well as the Bethe Ansatz results for the superalgebra Osp(n\2m). 
We remark that much of our motivation, concerning such general picture, was prompted by 
our previous effort of presenting the algebraic Bethe Ansatz solution [[HJ for the 4-dimensional 
representation of the supersymmetric spl(2\l) vertex model 0, [17], [l^. 19]. 



We have organized this papers as follows. In the next section, in order to make this paper 
self-contained, we review the basic properties of the braid-monoid algebra and its rational 
Baxterization. A convenient representation for Bethe Ansatz analysis is then presented for 
the Osp(n\2m) symmetry. In section 3 we formulate the eigenvalue problem for the transfer 
matrix in terms of the quantum inverse scattering method. The fundamental commutation 
rules are explicitly exhibited. In section 4 we elaborate on the construction of the eigenvectors 
and eigenvalues of the B n , C n , D n , Osp(2n - 1|2), Osp{2\2n-2), Osp{2n-2\2) and Osp{l\2n) 
vertex models and in section 5 the corresponding nested Bethe Ansatz equations are derived. 
These results allow us to conjecture the Bethe Ansatz equations of the general Osp(n\2m) 
chain. Section 6 is reserved for our conclusion and remarks on the universal picture we have 
found for braid-monoid vertex systems. In appendices A and B we present details concerning 
the two and the three particle state, respectively. In appendix C we collect some useful relations 
for the supersymmetric formulation of the quantum inverse scattering method. 



1 We remark that the algebraic Bethe Ansatz solution for the D n vertex model has been previously discussed 
by de Vega and Karowiski in rcf. []l4| by using a different construction than the one pursued here (see further 
discussion in section 5) . 
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2 The rational braid-monoid solution 



The braid-monoid algebra [[!]] is generated by the identity /, the braid operator bi and the 
monoid operator Ei acting on the sites % of a chain of length L. In general, these operators 
satisfy a set of relations which goes by the name of Birman-Wenzel- Murakami algebra |20], |J. 



We recall that for a braid operator 6, we mean a object satisfying the braiding relation ( 
bibi±\bi = bi±\bibi±i ) such that b{bj = bjbi for \i — j\ > 2. In this paper, we are interested in a 
degenerated point of this algebra, when the braid operator bi and its inverse b~ x are identical, 
i.e. bi = bj~ . Here we choose the braid operator as the graded permutation operator bi = Pf, 
defined by the following matrix elements [|ll| 

[p?)i = (-i) p(oM6) « c (i) 

where p(a) is the Grassmann parity of the a-th degree of freedom, assuming values p(a) = 0, 1. 



This is a generalization of the standard operation of permutation |TT[ , which distinguishes the 



'bosonic' (p(a) = 0) and the 'fermionic' (p(a) = 1) degrees of freedom. The monoid Ei is a 



Temperley-Lieb operator [£T| and satisfies the relations 

Ef = JCE t 

EiEi±iEi = Ei (2) 

E i E j = E j E i ; \i-j\>2 
where /C is a c-number. The choice we have made for the braid operator ( see equation (1) 
) greatly simplifies the constraints between the braid and the monoid. It is possible to show 
that the constraints closing the braid-monoid algebra on such degenerated point are given by 
( see e.g refs. [E 



PfEi = EiPf = tEi 

, (3) 

EiPl±\Pi = Pi± 1 PiE i ± 1 = EiE i±1 
where the constant i assumes only the values ±1. Any other constraint coming from the 
Birman-Wenzel-Murakami algebra [I], can be derived from (3), from the braiding properties 
of P 9 , the fact that {Pf) 2 = I, and the Temperley-Lieb relations (2). Lastly, we note that 
these set of relations are invariant by the transformation t — *■ — t and Pf — > —Pf . 



It turns out that the algebraic relations (2) and (3) can be 'Baxterized' in terms of rational 
functions. In other words, it is possible to find a solution R(X) of the Yang-Baxter equation 
in terms of certain combination of the identity, Pf and Ei. The solution comes in terms of 
rational functions [^, |3|] and is given byQ 

A 



R(X) = I + XP 9 - 



X- A' 



-E 



(4) 



where A = . This solution has a quasi-classical analog^, it is regular at A = and 

also has a crossing point at A = A. The next step is to search for explicit representations of 
the Temperley-Lieb operator Ei. In order to satisfy the braid-monoid restrictions (2,3), one 
possible choice is to set the following Ansatz for the monoid Ei 



(5) 



where a a b are the elements of an invertible matrix. A quite general representation for matrix 
a can be found in the context of the Osp(n\2m) symmetry. The integer n and 2m stands 
for the number of bosonic and fermionic degrees of freedom, respectively. The superalgebra 
Osp(n\2m) combines the orthogonal 0{n) and the sympletic Sp(2m) Lie algebras ( see e.g. 
ref. GHl ), and its element Z satisfies 



Z + M 0sp Z st Mo l sp = (6) 
where the symbol st denotes the supertranspose operation and the matrix Mo sp is given by 



M, 



Osp 



O 



o 



2mxn 



V 



/ 



V 



nx2m 



O m y.m Im x m 



\ 



(7) 



In this paper the i?-matrix is read as -R(A) = i?(A)^e ac ® ey, where the matrix elements of e a b are 



abed 



3 The quasi-classical r-matrix can be obtained by redefining A — ► — and by expanding it around r\ = 0. 
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where Ikxk and Okxk are the identity and the null k x k matrices, respctively. The elements of 
matrix (7) can be used as an explicit representation |22j for a a b, i-e. a a b = [Mo sp \ab, and the 
Temperly-Lieb parameter /C is then fixed by 

K = n - 2m (8) 

For general values of n and m, however, such representation breaks the U(l) invariance of the 
monoid [24] [], and it is not appropriate for Bethe Ansatz analysis. Usually, the lack of U(l) 



invariance induces extra difficulties on the formulation of a Bethe Ansatz, and therefore we 
should look for other alternatives. This problem can be resolved by the following construction. 
The monoid preserving the symmetry U(l) is built in terms of an anti-diagonal matrix a, 
whose elements are either +1 or —1. The integer m is the number of minus signs (—1) and n is 
the anti-trace of a. Furthermore, concerning the grading structure, the elements a a b only link 
degrees of freedom of the same specie, namely p(a) = p(b). The possible ways to distribute ±1 
in the anti-diagonal, for a fixed K, = n — 2m, are then related by permutation of the grading 
indices (canonical transformations). However, one needs to make sure that t, the elements 
a a b and the parities p(a) satisfy the braid-monoid relation (3), that is ia a b = (— l) p ^a&a- For 
example, one possible representation for matrix a is the following block anti-diagonal structure 



&Osp(n\2m) 



' O O T ^ 

^mxm ^mxm -^nxn 

O m xm -^m X m O n xn 

\ 1-mxm O m xm O n xn j 



(9) 



where Ikxk is a k x k anti-diagonal matrix. In this case, the two compatible sequences of 
grading are f x ■ ■ ■ f m h ■ ■ ■ b n f m+1 ■ ■ ■ f 2m for i = 1 and b 1 --- b m j x ■ ■ ■ f n b m +i ■■■b 2m for i = -1. 
One extra advantage of our construction is that the vertex models B n , D n and C m can be 
nicely represented in terms of the limits m — > and n — > 0, respectively. More precisely, we 
have 

I O m xm T-m / !/■ : - 
&B n = 22n+lx2n+l) a D n = -^2nx2n, «C m = | | (10) 

~^m,xm, Om.xrr 



4 We remark that the Osp(l\2m) is an exception to this rule. In this case, a canonical transformation can 
bring (7) in a U{\) invariant form. 



In this paper, besides the B n , C n and D n models, we are primarily interested in the super- 
symmetric models related to the Osp(2n — 1|2), Osp(2\2n — 2),Osp(2n — 2|2) and Osp(l\2n) 
symmetries. The first reason is because, as we shall see below, the nested Bethe Ansatz for- 
mulation for the Osp(2n — 1|2), Osp(2\2n — 2) and Osp(2n — 2|2) models goes fairly parallel 
to that of the B n , C n and D n vertex models, respectively. Secondly, because they exhaust the 
main basic symmetries present in the general Osp(n\2m) superalgebra. In order to see these 
relations, we show in figure 1 the Dynkin diagrams of the superalgebra Osp{n\2m) as well as 
those of the Lie algebras B n , C n and D n . We notice that the Osp(l\2n) superalgebra has a 
special structure of roots. In fact, such special character will be present in many points of our 
Bethe Ansatz analysis of the Osp{l\2n) vertex models. 

We turn now to the analysis of the Boltzmann weights of these vertex models. In general, 
these are multistate vertex models having one of q possible states on each bond of the two 
dimensional square lattice. The functional form of the Boltzmann weights depends directly 
on the values of i, /C. For the Osp(n\2m) vertex models, the weights also depend on the 
sequence of grading that has been chosen. In table 1 we have collected the values of q, i 
and K for the vertex models B n , C n , D n , Osp(2n - 1|2), Osp{2\2n - 2), Osp(2n - 2|2) and 
Osp{l\2n). Moreover, for the first three supersymmetric models we have used the grading 
fb- ■ -bf and for Osp(l\2n) we have taken the sequence b ■ ■ ■ bfb ■ ■ ■ b. The basic reason for 
this choice is because it uses the minimum number of fermionic degrees of freedom in the 
grading sequence compatible with number /C and our choice of t (see table 1), and therefore 
saving the commutation rules of the presence of many extra minus signs. We remark that the 
structure of equation (9) is compatible with the above mentioned grading sequences for the 
Osp(2n — 1|2), Osp(2n — 2|2) and Osp(l\2n) vertex models. For the Osp(2\2n — 2) system, 
however, we still have to perform an extra transformation in matrix (9), by exchanging the least 
— 1 and the (n + l)-th +1 degrees of freedom. More specifically, we have that cto sp (2\2n-2) = 
anti — diagonal(li, • • • 1„, — l n +i, — ln+2, • • • , —^2n-i, ^2n)- This transformation is quite helpful 
because it allows us to treat the nesting problem for the C n and Osp(2\2n — 2) models in a 
common way. 
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Taking into account these considerations, we find that there are only few possible distinct 
functional forms for the Boltzmann weights. In fact, as we shall see in next section, our 
algebraic formulation will require at most five distinct functional dependences. We name these 
Boltzmann weights by a(A), b(X), c n (X), d n (X), e n (X) and they are summarized on Table 2. To 
some extent, this is the miraculous fact of integrability, which in our approach is encoded in 
the commutation rules to be presented in next section. 



3 The eigenvalue problem and the commutation rules 

We start this section by describing the eigenvalue problem for the transfer matrix T(A) associ- 
ated to the vertex models defined in section 2. We are interested to determine the eigenvalues 
and the eigenvectors of T(X) on a square lattice of size L x L. The diagonalization problem is 
defined by 

T(A) |$) = A(A) |$) (11) 



An important object in the quantum inverse scattering method |7|, §, [TO], |TTJ is the Yang- 



Baxter algebra of the monodromy matrices. The monodromy matrix T(A) acts on the tensor 
product of an auxiliary space A = C q and on the quantum Hilbert space C Lq . The transfer 
matrix T(A) is the trace of the monodromy matrix T(A) over the auxiliary space A, i.e., 
T(A) = Tr^T(A). A convenient way of writing the monodromy matrix is in terms of a product 
of operators Cm(X) as 

T{X) = £al(X)£al-i(X)....£ a1 (X) (12) 

where £^j(A) are q x q matrices acting on the lattice sites i = 1, • • • , L whose elements are 
operators on the Hilbert space C Lq . The elements of the vertex operator £^(A) are related to 
those of the i?-matrix (4) by a permutation on the C q x C q tensor space as 

£&M = fl£(A) (is) 

and the Yang-Baxter algebra for two monodromy matrices with distinct spectral parameters 
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reads 



R(X - fj)T(X) T(fj) = T{n) ® T{\)R{\ - n) (14) 

The intertwining equation (14) is the corner-stone of the quantum inverse scattering ap- 
proach, allowing us to derive the fundamental commutation relations. Unfortunately, there 
is no recipe to immediately find the appropriate commutation rules from the Yang-Baxter 
algebra. Much of the insight comes from the properties of the vertex operator £^i(A) itself 
and from the reference state we choose to begin the construction of the Hilbert space. As a 
reference state we take the standard ferromagnetic pseudovacuum given by 

/1\ 



io> = n®i°>i. i°>i 







(15) 



where q stands for the length of the vectors |0).. The vertex operator £^j(A), when acting on 
state 1 0)^, has the following triangular property, i.e., 

( a(A) |0) i * * * * * \ 

&(A) |0>. ... * 
6(A) 10)^ ... * 



Aai(A) |0> 4 = 



V 























(16) 



.. 6(A) 10), 

e n (X)\0)J qxq 

where the symbol * represents some values that are not necessary to evaluate explicitly for 
further discussion. The next step is to write an appropriate Ansatz for the matrix represen- 
tation of T(A) on the auxiliary space A. The triangular property (16) suggests us to seek for 
the following structure 



T(A) = 



/B(X) B(X) F(X)\ 

(7(A) i(A) B*(X) 
\C(X) C*(X) D(X)J 



(17) 



qxq 
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— * — » — * — * 

where B(X), B*(X), C(A), C*(A) are two component vectors with dimensions 1 x (q — 2), 
(q — 2) X 1, (g — 2) x 1, 1 x (g — 2), respectively. The operator A(A) is a (g — 2) x (g — 2) matrix 
and we denote its elements by A a b(X). The other remaining operators B(X), C(X), F(X), and 
D(X) are scalars. Putting them all together, we then have a (q x g) matrix representation 
for the monodromy matrix T(X). Taking into account these definitions, the diagonalization 
problem (11) for the transfer matrix T(A) becomes 

[B(X) + £ A aa (A) + D(A)] |$) = A (A) |$) (18) 

a=l 

As a consequence of definitions (12,17) and the triangular property (16), we find that the 
diagonal operators of T(X) satisfy the following relations 

B(X) |0> = [a(X)] L |0> ; D(X) |0> = [e n (A)] L |0> ; A aa (X) |0) = [b(X)] L |0> , a = 1, . . . , q-2 (19) 

as well as the annihilation properties 

C a (A)|0) = 0; C a *(A)|0)=0; C(A) |0) = 0; A ab (X) |0) = 0, (a, 6 = 1, ■ • • , q - 2; a ^ 6) (20) 

This means that the reference state is an exact and trivial eigenvector with eigenvalue 

A(A) = [a(X)] L + (g - 2)[b(X)} L + [e n (X)} L (21) 

and also that the fields -B(A), B*(X) and F(X) should play the role of the creation operators 
on the reference state. In order to construct the full Hilbert space we need to find the com- 
mutation relations between the creation, diagonal and annihilation fields. In principle, all the 
information concerning commutation rules are encoded in the integrability condition (14). The 
basic problem is to collect them in a convenient form. For instance, for the 6-vertex model 
[0, H] and its multi-state generalizations [JT0|, [1^, [13| they come almost directly, after substitut- 



ing the appropriate form for T(A) and the associated i?-matrix on the Yang-Baxter algebra 
(14). For the models we intend to analyse in this paper, that is the rational braid-monoid 
vertex models (4), this is not the case and some additional work is necessary. For example, in 
order to get the 'nice' commutation rule between the creation operator B(X) and the diagonal 



field A(X) we also have to disentangle the commutation relations between B(/j,) and B*(X). 
The trick goes much along the lines we have already explained in details for the spl(2\l) su- 
persymmetric vertex model ||15|| . It turns out that these ideas can be generalized for the B n , 
C n , D n , Osp(2n — 1|2), Osp(2\2n — 2), Osp(2n — 2|2) and Osp(l\2n) multi-state vertex models, 
and therefore we have been able to derive their fundamental commutation rules. As stressed 
below, the main procedure is quite cumbersome, and here we just list our final results for the 
most important commutation rules. Between the diagonal operators A(X), B(X), D(X) and 
the creation field B (A) we have 

1 r -, , ~m„ , 1 



A(X) <g> B(fi) 
d n (X - fi) 



■[B0.)«A(A)].*«(A-rt- 



B(\)®A(ji) + 



e n (X - n) 



B*(X)B(») + _J_F(A)d(Ai) - {1 Jf (A / )} F(/i)C(A) 



6(A- M ) 

a(/x — A) 
= % - A) 



5(/x)B(A) 



6(A-ai) 
1 



6(A* - A) 



(22) 
(23) 



D(X)B(fi) = -^j- — ^-B(fi)D(X) + 
e„(A - /ij 

c n (A - /i) 



1 



e„(A - fi) 



F(fz)C*(X) 



(24) 



e„(A-yu) ' e n (A-/i) 

where X^ 1 -* (A) is a factorizable auxiliary i?-matrix responsible for the first nested Bethe Ansatz 
structure. In what follows we shall see that it is also useful to introduce a second factorizable 
i?-matrix JW(A). The matrix elements of and JW(A), however, are related to each 
other under permutation of their horizontal and vertical spaces, namely XW(A)^ = X^(X)^. 
This distinction is necessary in order to include the Osp(l\2n) model in our discussion. This 
model is an exception because its auxiliary i?-matrix X^'(X) is no longer T-invariant. All the 
other vertex models, however, have the property A^'(A)^ = X^'(A)^, and therefore X^\X) 
and X«(A) are indeed identical. In table 3 we have collected the Boltzmann weights of the 
iZ-matrix X^(A) on a certain k-ih step of the nested Bethe Ansatz. We also show the value 
of crossing parameter corresponding to X^ k \X). From this table we notice that the pairs 
of models: {B n ,Osp(2n - 1|2)}, {C n , Osp(2\2n - 2)}, and {D n , Osp(2n - 2|2)} share the 
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same auxiliary X( k \X) matrix. In fact, this is not that surprising if one takes into account 
the similarities between their root structure ( see figure 1). Furthermore, for these models, the 
matrix is precisely the i?-matrix associated to the -B„_fc, C n ~k and -D„_fc vertex models, 
respectively. However, since the vertex models on a given pair have distinct Boltzmann weights, 
we shall see that their eigenvectors and eigenvalues will be in fact different too. As before, 
the Osp(l\2n) model is an exception. Here, the nesting still keeps the same structure as the 
original model and we have the embedding Osp(l\2n) C Osp(l\2n — 2) C •■■ C Osp(l\2). 
Finally, the vector £ and its "conjugated" £* are given in terms of the matrix a, defining the 
monoid operator present in the i?-matrix They are defined by 

q-2 q-2 

I* = E a ab(e a ® e b ), = ]T a ab(e a ® e b ) (25) 

a, 6=1 o,6=l 

where denotes the elementary projection on the i-th position. 

We now will give other important commutation relations. The commutation rules between 
the scalar creation operator F(X) and the diagonal fields are given by 

' r 5(A) ®B*(/j) -B*(X) ®B(ii)} (26) 



b(X - fj) 



B(\)Ffr) = ^-^-F^)B(X) - *!<JLVf(\)BM - ^m?.{5(A) ® B(»)} (27) 
e„(/i-A) e n (fi-X) e„(/i-AJ 

D(X)F^) = ^p±F(»)D(X) - 5£^F(\)DM - ^^MB*(X) ® £*(/,)} (28) 
e n (A-/i) e n (X-n) e n (X - /i) 

and the commutation relations between the creation operators are 

B(X) ® B(jm) = 1 [gfr) ® B(A)].lW(A - /i) 
a(A — /i) 



d n (A - /i) 
e n (A - //) 



^)B(,)-H^0} WA) 



(29) 



P(»).F(P)] = (30) 
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fl(A)F(/i) = — i_B(/i)F(A) + *^-^F(a*)J3(A) (32) 

We close this section with the following remark. In this section, we have kept our presen- 
tation concerning the basic properties of the quantum inverse scattering approach as general 
as possible. All the vertex models have been treated by the same and standard way of formu- 
lating the quantum inverse scattering method. We notice, however, that the solution of the 
super symmetric Osp(n\2m) models could also be presented in terms of a graded framework 
from the very beginning. This is possible because the i?-matrix (4) has a null Grassmann 
parity, and consequently can produce a vertex operator Cjn satisfying either the Yang-Baxter 
equation (standard formulation ) or its graded version [[11], [J. The last choice, however, is 
the most natural way of formulating the vertex operator C^i for the Osp(n\2m) models, if 
one wants to make a real distinction between the bosonic and fermionic degrees of freedom. 
In other words, the graded quantum inverse method makes sure that the fermionic degrees 
of freedom anticommutes even if they act on different lattice sites. In order to accomplish 
this "non-local" property, one has to use the supersymmetric (graded) formalism developed in 
refs. [j5], 1 1 ]| ( see also refs. [29| ), which basically changes the trace and the tensor product 
properties by their analogs on the graded space. In appendix C we have summarized this ap- 
proach for the Osp(2n — 1|2), Osp(2\2n — 2), Osp(2n — 2|2) and Osp(l\2n) vertex models. We 
notice, however, that the supersymmetric formulation does not simplify the original problem 
of the diagonalization of the corresponding transfer matrices. We shall come back to this point 
again in section 5, where the final results for the eigenvalues and Bethe Ansatz equations are 
presented. 



4 The eigenvectors and the eigenvalue construction 

In the quantum inverse scattering scheme the eigenvectors are constructed by acting the cre- 
ation operators on the ferromagnetic pseudovacuum |0). The excitations over the pseudovac- 
uum |0) present a multi-particle feature, characterized by a set of variables {A x , ...,A^} 
which are determined after a posteriori analysis. The way we are going to start our discussion 
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is much inspired from our previous results for the supersymmetric spl(2\l) vertex model |15[ . 
Indeed, the construction we shall present here is a non-trivial generalization of ideas discussed 
by us in ref. ||15| . Due to the presence of many kinds of creation fields, it is convenient to 
represent the mi-particle state ^ mi (X^\ • • ■ , X$)) by the following linear combination 



$ mi (Al 1) ,---,AS)) = $ mi (AS 1) ,---,AS)./|0) 



(33) 



where $ mi (A^ ,---,A^) and f are multi-component vectors with (q — 2) mi components. 

(i) 



The dependence on the creation fields is encoded in the vector $ mi (A] 1 '', ■ ■ • , \$.) and T 



is a 



constant in this space of fields. We shall denote the components of T by !F am i- - 01 . 

The simplest excitation, i.e., the one-particle state, can be built only in terms of the creation 
fields B(X), namely 

$ 1 (\?)=B(\W) (34) 
and as a consequence of equations (33,34), the one-particle state is given by 



$i(A 



B a (X?)F a \0) 



(35) 



where here and in the following repeated indices denote the sum over a = 1, ■ • • , (q — 2). The 
commutation relations (22 — 24) and properties (19, 20) can be used to solve the eigenvalue 
problem for ^(A^)). The solution is based on the following relations 



b(x) $ 1 (aS 1) )) = ^)!^ka^ $ i( aS 1} ) 



6(aS 1} - A) 



1 



a(AS 1) )] L J B a (A)F a |0) (36) 



K\? ~ A) 



D(X) ^(X?) 



en(A-A^) 



9-2 



$> aa (A) $ x (Ai 



a=l 



b(X-X?) 



e n (X-X?) 



l + ( g -2)(A-A? ) )- 



^HX\ L, )] L UB:(X)F b \0) (37) 



(A - A?>) 



[b(X?)] L B a (X)F a \0) + 



dn{X-X^) ^ n{1) 

e„(A - X?) 



(A- aS 1] + A«) 
«(Ai i; )] L a5:(A)F"|0) 



[b(X)] L ^(X? 

(2 



b(x - aS 1} ) 

The terms proportional to the eigenstate $i(A ( j 1 ' ) )^ are denominated 'wanted' terms and 
contribute to the eigenvalue A(A, A^ ). The remaining ones are the so called 'unwanted' terms 
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and they have to be canceled by imposing further restriction on variable A^ . Such constraint 
goes by the name of Bethe Ansatz equation. From relations (36 — 38), for the one-particle 
state, we find a single equation given byQ 

L 

= 1 (39) 



a(X\ 



(ih 



MAS"; 

The feature of the two-particle state is a bit more complicated. Now, the scalar operator 
F(X) starts to play an important role on the eigenvector construction. This becomes clear if 
one takes into account the commutation rule between two creation fields of type -B(A) (see 
equation (29)). For the two-particle state we have to seek for a combination between two fields 
of type -B(A) with the scalar field F(X). The structure of commutation rule (29) suggests us 
to take the following combination for the vector ^(A^, A^) 

$ 2 (aS 1} , A?>) = B(X?) ® S(A«) + h(X?, xP)F(X?)B(X l 2 )g (40) 

One way of fixing the function h(Xi , A^) is to notice that the eigenvalue problem for the 
two-particle state generates certain kind of 'unwanted' terms which are proportional to 
They are given by 

F(X)D(X?)B(X?); B(X).B*(X?)B(X?); f.[if*(A) ® &(\P)]B{)$)) (41) 

We call these structures 'easy unwanted' terms, because they can be automatically canceled 
out by an appropriate choice of the form of the function ^(A^, A^), namely 

fc(A<", A<") = MA« - A? ») = (42) 

e n (Ai - Xy) 

There is a more elegant way, however, to determine the function /i(A^ 1 ' ) , A^). In general, 

we expect that the extra constraints on variables A^ , i.e. the Bethe Ansatz equations, are 

invariant under index permutation. This means, for the two-particle, that vector $ 2 (^i \ A2 ) 

and ^(A^, Ai ) have to be related in some sense. The commutation rule (29) itself gives us a 
5 This means that the degeneracy of the one-particle state is (q — 2), because no further constraint is 
necessary for vector J- ' . 
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hint how these two vectors can be related to each other. In fact, it is not difficult to conclude 
that the following exchange property 

* 2 (A« A«) = * 2 (A« AS") X 'yj' -£> (43) 

is valid, provided the function ^(A^, A^) is fixed as in equation (42). To prove (43) we have 
used the remarkable identity 

fc(A« A<«) 

which relates vector £, the auxiliary matrix and the Boltzmann weights. This is another 
way to see the role of vector £; it is an eigenvector of the auxiliary matrix with defined 
eigenvalue where the function h(X^\ A^) enters in a symmetrical way. 

In order to cancel out other remaining 'unwanted terms', it is necessary to impose further 
restriction on variables {A^, A^}. In appendix A we show how many different kinds of 
'unwanted terms' can be canceled out by the following Bethe Ansatz equations 



eX«(A« - A?>) = ^'^ {1 - tb(X? - A?>)}< ( I I) 



a(X\ 



0-h 



(i)> 



L 



II rf-Af'j^^r^Ptf-Aj 1 ')^, i = l,2. (45) 



" " * J > w,(l) ,(l)x 



The same kind of reasoning can be used in order to construct the 3-particle state. Here, 
however, it is already interesting to write it in terms of certain recurrence structure. Inspired in 
our previous construction for the supersymmetric spl(2\l) model, we begin with the following 
Ansatz 



*3(A«,A«,A«) = 5(A«)0$ 2 (A«,A«) + 



F(X?)B(X?)B(X?)\ hW, Xf>, + [£® F(XY>)B(X^)B(X^)\ h 2 (X\\X^, X^>) 

(46) 

As before, the functions hi(Xi~\ X^\ X<p) and h 2 (Xi~\ A^\ A^) can be fixed either by col- 
lecting together the 'easy unwanted' terms or by using the permutation symmetries A^ — > A^ 



>(') \W \W 



a(1) ,(lh 
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and Ai — > Xp. We have found that they are given by 

h f\W X 1 X ^ - d "l A l - A 2 ) Q l A 3 - A 2 j j (A7 , 

hl{Xl 2 3 "^-a^ka^-aW) 7 (47) 

" 2(Al 2 "^-A^KA^-Af^ 23 ( 2 ~~ 3 } (48) 

While the permutation A 2 — *• A3 is easily verified, that concerning the symmetry Aj — > 
Ag is quite cumbersome. One of the main difficulty, for example, is that we have also to 
commute B(X^) and F(A 2 )■ Therefore, the commutation relations (31 — 32) between the 
creation fields play an important role to disentangle the permutation Aj 1 — > A 2 . Besides that, 
certain additional properties between auxiliary matrix X^ 1 ', vector £ and the field B(X) are 
also necessary For sake of completeness, we have collected the details of the analysis of the 
three-particle state in Appendix B. Considering the structures of the two and three-particle 
state (equations (40) and (46)), it becomes clear that the mi-particle state can be generated 
by means of a recurrence relation. The general mi-particle state can be obtained by using an 
induction procedure, and our final result is given byf] 

n(\ {1) )(y>$ \M) V n v 1 ~ A i j TT a l A fc ~ A j j 

j=2 e n {A 1 - Aj j k=2,kytj 0{A k - Aj J 



f ® F(A«)4- 2 (A«, • • • , A«, A«, . . . , AW 1 )5(A«)] n ^±^— ^ 



l 2 "(A^-Af) 



(49) 



and under a consecutive permutation A^ — > A^ the mi-particle state satisfies the following 



relation 

j,j+n A j in 
a( A f " A fft) 



<B A (1) A (1) A«^ - $ fA (1) A (1) A (1) \(D\ ^id+n A j /r m 
^miV A l > ) A j ) A i+1) ) A mi) — v milA > > A/+l> Aj , •••,A V 7 /J j^-^ ^jy^ l du J 



3 The indices under Xj, l + i(X k Aj ) indicate the positions on the Hilbert space where this matrix acts in 
a non-trivial way. 
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The following remark is now in order. So far we have not commented about the role of the 
creation field B*(X). From the integrability condition (14), we can verify that the commutation 
rules between the creation field B* (A) with the diagonal operators and with the scalar creation 
field F(X) are similar to those in equations (22 — 24, 26 — 32). Basically, we have to change B(\) 
by B*(X), to replace the definitions of A(X) and £ by their transpose ones, and to interchange 
the diagonals operators B(X) <-> D(X). The three creation fields do not mix all together, but 
only in pair as {B*(X),F(X)} and {B(X), F(X)}. We shall then call B*(X) the " dual " of B(X). 
This dual creation field together with F(X) can also be used to built eigenvectors, and by using 
the replacements mentioned above we obtain an expression analogous to that of equation (49). 
We think that these two possible ways of constructing the eigenvectors go back to the fact 
that the vertex operator Cj^{X) becomes naturally crossing invariant after multiplying it by 
the factor (A — A). Indeed, it is not difficult to see that their corresponding eigenvalues are 
related to each other by performing the crossing shift x — > A — x both in the transfer matrix 
parameter and in all the Bethe ansatz rapidities. It seems interesting to analyse to what extent 
we can benefit from this property of the eigenvectors in order to provide new physical insights 
for these vertex models. 

Let us now turn to the Bethe Ansatz conditions on the variables {A^}. A typical unwanted 
term, coming from the eigenvalue problem (18), arises when the variables X^p of the mi-particle 
state is exchanged with the transfer matrix rapidity A. For instance, this is the case of an 
unwanted term having the structure B ai (A) B a2 (X l 2 '')••• B a (A^). This unwanted term is 
produced by the action of the diagonal operators B(X) and J2 a Am (A) on the mi-particle state. 
From the commutation rules (22, 23) we are able to show that they can be collected in the 
following forms 

- ^P 1 ^ 5 S'S ^ 1 '^'^^' ' " B "»>< A 2> l°> (51 > 

ft (A ~ X ) 3=2 0{\j - \ ) 

and 

fl kn m . 01 , r"'(A = Ai 1 ', {^})Z:Z;^ --^.WiWA?') ■ ■ -B^l) |o) 

°l A l - A J 3=2 0(A - Aj ) 

(52) 
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where T^(X, {A^ }) is the transfer matrix of a inhomogeneous vertex model defined in terms 
of the auxiliary .R-matrix X^'(X) by 

t (i) (a, = - x^wai^ix - x^tz ■ ■ ■ [x«(a - x^txr (53) 

Similarly, the same kind of reasoning can be pursued for any A$ , thanks to the prop- 
erty (50) which relates vector $ mi (Ai , • • ■ , A^) to $ mi (A^ , ■ • • , XQ) by cyclic permutations. 
Thus, from equations (51) and (52) the unwanted term B ai (X^) . . . B ai (X) ■ ■ --Ba mi (A^) can 
be canceled out provided that 

t«(a = a« {xf})Z; a b :ir b ^-- b \i = i, . . . , mi (54) 



a(X 



This generalizes equation (45) for an arbitrary value of the number of 'particles' rri\. In 
general, the action of the diagonal operators on the mi-particle state generates many other 
species of unwanted terms. A systematic way to collect all of them in families of unwanted 
terms with a defined structure for a general value of mi has eluded us so far. However, we 
remark that for those we have been able to catalog, such as B^X^) . . . B a .(X) ■ ■ ■ B am (A^), 
condition (54) has been explicitly verified. Furthermore, the explicit checks we have performed 
in two and three-particle states leaves no doubt that restriction (54) is the unique condition]] 
to be imposed on the variables {A^, • ■ ■ , A^}. Anyhow, a rigorous proof should be welcome, 
because probably it will shed extra light to the mathematical structure we have found for the 
eigenvectors. 

By the same token, the eigenvalue A(A, {A^ }) of the mi-particle state can be calculated 
by keeping only terms proportional to the eigenstate $ mi (A^, • • • , A^)V For instance, that 
proportional to the term B ai (X^) . . . £? ai (Aj- ) • • ■ B^ (A^) can be determined by the commu- 
tation relations (22 — 24). By keeping the first terms each time we turn the diagonal operators 
7 In a factorizable theory, which is our case here, the two-particle structure already contains the main flavour 
about the Bethe Ansatz equations ( see equation (45) ). 
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— * 

over one of the B(X) components, we find that 

mi i 

^] L n 77^(1)7 AW (A, {AS 1} }) (55) 
i=i 6(A- A^ 

where AW(A, {A- 1 ^}) is the eigenvalue of the inhomogeneous transfer-matrix T^(X, {A- 1 ^}), 

T (1) (A, {aS^})^ 1 .:::^ 6 ---- 61 = A (1) (A, {A^}).^— (56) 

and therefore the Bethe Ansatz equations (54) are then disentangled in terms of the eigenvalue 
A«(A,{Al 1) }) as 

H - Af) ^-^ = A«(A = xf\ {Af }), < = l > ...,m 1 (57) 

This completes the first step of our analysis, because we still need to find the eigenvalue 
A (1) (A, {Aj 1} }) of the transfer-matrix TW(A, {Af ) }). This is the so-called nested Bethe Ansatz 
problem and we shall deal with it in next section. 



a(A 



ah 



5 The nested Bethe Ansatz problem 



In the last section we were left with the problem of diagonalization of the inhomogeneous 
transfer matrix T^(X, {A^}). This problem can still be solved by the quantum inverse scat- 
tering approach, once the Yang-Baxter algebra (14) is extended to accommodate the presence 
of inhomogeneities. Formally, the corresponding monodromy matrix can be written as 



t«(a, {a; i} }) = 6Zm - a2)4Vi(a - a£!-i) • • -4i(A - aS : 



0-h 



(58) 



where the components of the vertex operator [£^(A)]^ are related to that of the i?-matrix 
[XW(A)]^ by the standard permutation 



cd 



£ (1) (A) = X (1) (A) 

ab 



cd 
ba 



(59) 
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In order to go on, we need first to analyse the properties of the auxiliary matrix jW(A). 
For instance, we have to verify if the triangular property of C^(X) still has the same structure 
as that present on the original vertex operator C(X) we started with ( see equation (16) ). In 
other words, we have to check whether or not the basic structure of the Ansatz (17) for the 
monodromy matrix is still appropriate. As we shall see, this depends much on the type of the 
original vertex model we are diagonalizing and also on their number of states per link q. We 
recall that in the first step we already lost two degrees of freedom, remaining (q — 2) states per 
link to represent the Boltzmann weights of matrix X^(X). Hence, if (q — 2) < 3 the Ansatz 
(17) certainly is not more the convenient one. 

Let us, for the moment, suppose that the i?-matrix X^(X) keeps the basic properties of 
the original model. In this case, we have basically to adapt our main results of sections 3 and 
4 to include the inhomogeneities {A^, . . . , A^}. For instance, the pseudovacuum CK 1 ^ is the 
usual ferromagnetic state, but now with length (q — 2) 



mi 



( i )) = o 



(i) 



i=i 







(60) 



q-2 



and the vertex operator Cj^(X — A- iJ ) has the following triangular property 



^a(A-A< 1) )|0< 1 )) i 

b(\-x ( i 1) )\oW). 



^(A-A« 1) )|0«) i 



6(A- Af>) |oW). 

^(A-A^joW) / ( ,_ 2)x( ,_ 2) 

(61) 



As before, if we assume the following structure for the monodromy matrix T^(X, {A^}) 
(BW(X,{XV}) BM(X,{X?}) f«(a,k (1) })\ 

(62) 



TU(X,{XM}) 



CM(X,{XV}) AW(A,{A«}) 5* (1) (A,{A, (1) }) 



VC«(A,{A«}) C^\X,{\? } }) ^ (1) (A,{A, w })y ( ^ 2)x( ^ 2) 



(!)• 
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we find that the operators of the monodromy matrix T^(X, {A- 1 ^}) satisfy the following 'di- 
agonal' properties 



mi 



mi 



3«(A, {A«}) |o«) = n«(A- A 4 (1) ) |0«) ; D^(X, {A«}) |o«) = e»-i(A - A«) |o«) ; 

i=l i=l 

mi 

42(A, {AS 1} }) |0«) = I] 6(A - A, (1) ) |o«) , a = 1, . . . , q - 4(63) 
as well as the annihilation properties 



i=i 



^(MA^}) 0W) = 0; C:«(A,{AW}) 0«) = 0; 



C^A,^}) 0«)=0; A«(A,{A«}) 0«) = 0, (a, 6 = 1,-.., ? -4; a ^6) 



(64) 

In order to diagonalize the transfer matrix T^(X, {A- 1 - 1 }) we have to introduce a new set of 
variables {xf\..., A^}, which are going to parametrize the multi-particle states of the nesting 
problem at step 1. Evidently, the structure of the commutation rules (22 — 24; 26 — 32) and the 
eigenvector construction of section 4 remains precisely the same. We basically have to change q 
to (q — 2) in the Boltzmann weights, a given operator 0(A) by its corresponding 0^(X, {X^}) 
and to replace the parameters {A^, . . . , A^} by {X^, . . . , A^} in the eigenvector expression 
(49). It turns out that the eigenvalue expression AW(A, {A^}) and the corresponding Bethe 
Ansatz equations will again depend on another auxiliary inhomogeneous vertex model, having 
now (q — 4) states per link. Of course, we can repeat this procedure until we reach a certain 
step / where the underlying auxiliary i?-matrix X® lost the basic features presented by the 
original i?-matrix we began with. In general, by using this "nested" procedure we can find 
a relation between the eigenvalues A«(A, {Xf }, . . . , {A^}) and A(' +1 )(A, {A^ +1) }, . . . , {\f ] }), 
on the steps / and I + 1, respectively. We basically have to " dress " our previous result (55) 
with inhomogeneities, and to consider the appropriate Boltzmann weights of the step / we are 
diagonalizing. This relation can be written as 



AW (A; {A? },-.., {Af} 



mi 



b(X 



A}' 



m ' +1 n(\( l+1 ^ — \\ mi m ' +1 

nxA-A^n -k^ -( +n^-,(A-Af)n , l „ 

j=i j=i b{X) - A) 3=1 J= i e (n _ (A - A} ) 



+ 
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mi rn l+ i 

n 6(a - a«) n 



"i KX - Xf +1) ) 



A (m) (A;{A^},{A5},...,{A«}) 



(65) 



This last equation has to be understood as a recurrence relation, beginning on step zero. In 
order to be consistent, we identify the zero step I = with the eigenvalue of the original 
transfer matrix T(A) we wish to diagonalize. Therefore, we are assuming the following iden- 
tifications A(°) = A, {A' " 1 } = and m = L. Analogously, the Bethe Ansatz restriction on 
the variables {A? } introduced to parametrize the Fock space of the inhomogeneous transfer 
matrix T®(\, {A^ }, . . . , {A^ }) on step I is given by 

TT \ j ~ A i ) TT i/aH+D ,(i+l)N a l A i ~ A 



n 



,-= =i 6(Af +1) - A? } ) 



n 6(A 



(I+l) _ 



i=l i^j &(A 

A (m) (A = A (m )){A (m )}) ... ){A (i) }); J = 1 ,..., m;+1 



(66) 



We now will particularize our discussion concerning the nesting structure for each pair of 
models {B n , Osp(2n - 1|2)}, {C n , Osp(2\2n - 2)}, {D n , Osp(2n - 2|2)} and for the Osp(l\2n) 
vertex models. For the pair {B n ,Osp(2n — 1|2)}, it is not difficult to check that the nesting 
structure developed above ( see equations (65,66) ) works in any step. In this case, the last step 
consists of an inhomogeneous Bi vertex model and the respective i?-matrix acts on C 3 x C 3 
tensor space. The vertex operator C Bl (X) has the following structure 



/o(A) 























^ 





6(A) 





1 























ei(A) 





di(X) 





d(A) 











1 





6(A) 























di(X) 





A(A) 





rfi(A) 























6(A) 
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ci(A) 





di(X) 





e x (A) 























1 





6(A) 





V o 























a(A)/ 



(67) 



where the Boltzmann weights a(A), &(A), ci(A), di(X) and ei(A) are listed on table 2 for Bi, 
and /i(A) = 6(A) +Ci(A). Due to the isomorphism B 1 = 0(3) ~ SU(2) k=2 , this vertex model is 
equivalent to the isotropic spin-1 XXX model ||28||. Therefore, the B\ model can be solved either 
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by adapting the known results for the spin-1 XXX model |28| to include inhomogeneities or by 
applying the general approach we have developed in sections 3 and 4. In the latter case, we 
remark that our construction for the eigenvalues, eigenvectors and commutation rules reduce 



to that proposed prviously by Tarasov in the context of the Izergin-Korepin vertex model [27 



The final result we have found for the eigenvalue A Bl (X, {Aj}, {Hj}) of the inhomogeneous Bi 
vertex model is 



A + l 



A B > (A, {A,}, = n(A - A, + 1) n ^ r 

j= i V-j ~ A 



+ 



i=i 



jj (a ^ (A-^ + l/2) (A-^-l) | ^ (A-AQ(A-A<-1/2) Q A-^ + l/2 



i=i 



3=1 



(A - Mi ) (A - H - 1/2) 



i=i 



(A - A, + 1/2) 



?= l A _ ^ _ 1/2 
(68) 



where {Ai,---,A„} are the inhomogeneities and variables {/ij} satisfy the following Bethe 
Ansatz equation 



^ H - Xi + 1 



n 



/ij - //jfc + 1/2 



j = l,...,m. 



(69) 



i=l 



fJ-j - Aj fe=1 Mi /ij -fi k -l/2' 

On the other hand, the situation for the pairs of models {C n , Osp(2\2n — 2)} and {D n , 
Osp(2n — 2|2)} is a little bit different. In these cases, we can proceed by using the recurrence 
relation (65) and the Bethe Ansatz condition (66) until we reach the steps I = n—1 and / = n— 2, 
respectively. For the C n and Osp(2\2n — 2) vertex models the nesting level I = n — 1 corresponds 
to the diagonalization of an inhomogeneous transfer matrix which Boltzmann weights have the 
6-vertex symmetry. Indeed, at this level we have to diagonalize the inhomogeneous C\ system 
whose vertex operator £ Cl (A) is given by 

/l 




£ c i(A) = (A + l) 







A 2 

A+2 A+2 

2 A 

A+2 A+2 



o\ 




1/ 



(70) 



The diagonalization of the 6-vertex models on a irregular lattice has appeared in many 



different context in the literature (see for instance refs. ]I2, [L3|, [14], [15], |2"5fl). By adapting these 



23 



results, in order to consider the particular structure of £ Cl (A), we find that the eigenvalue 

A c > (A, {A,}, is 

a c - (a, {a,}, = n<A - a, + 1) n ^¥ + n (A ",^",V" n A " " 3 ' + 2 



i=l 



=1 A*i ~ A 



i=l 



(A — Aj + 2) 



=i A - N 



where variables {/ij} satisfy the Bethe Ansatz condition 

/ij — Xi 



n 



m . . 9 

nf*j f*>k a 
~ 7. ~> 



l,...,m. 



(71) 



(72) 



= i - Aj + 2 fc=1 fe _^. fij - fi k + 2 

We now turn to the analysis of the last nesting level / = n — 2 for the models D n and 
Osp(2n — 2 1 2). The final stage for these systems consists of the diagonalization of the D 2 
vertex modelQ. It turns out, however, that the D 2 weights can be decomposed in terms of 
the tensor product of two 6- vertex models. More precisely, the vertex operator C D2 (X) can be 
written as 

C D2 {\) = cl~ ver ® C 

where the two 6- vertex structures are given by 



6— ver 



(73) 
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X + l) 



(74) 



Consequently, the eigenvalues of the model D 2 are given in terms of the product of the eigen- 
values of the two 6- vertex models defined in equation (74). In the presence of inhomogeneities 



{A 



l) 



A n } we find that these eigenvalues are given by 



A D > (A,{A,},{/4} 



m + \ 11+ 1 
n A_ft_l +n 



A -A, 



m_|- 



\-i4 + i 



j i A — fi ~- fJi A — Xi + 1 A — [l- 

n \ , , — n(A-A,+i)+n / j „- ii(a-a, 



x 



A -//• 



'3 1=1 



i=i 



A — //,• 



i=i 



(75) 



8 Naively, one would think that this model can still be reduced to the D\ vertex model. This is not the case, 
because the i?-matrix of the the D\ vertex model (A = in equation (4)) is no longer regular at A = 0. 
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where variables and {^j} parametrize the multi-particle states of the inhomogeneous 

models related to C®~ ver (\) and C®~ ver (\), respectively. They satisfy the following Bethe 
Ansatz equations 



n 



\ m± 
~ n 



i4 - ^ - 1 
(4 - + 1 ' 



j = l,...,m±. 



(76) 



i=1 A*j "l I - fc=1 kj Lj 

Finally, it remains to investigate the Osp(l\2n) vertex model. In this case, the nesting 
recurrence relations (65,66) are valid for any level. On the last step we have to deal with the 
inhomogeneous Osp(l\2) model, possessing the following vertex operator 
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where the Boltzmann weights a(A), 6(A), Ci(A), di(X) and ei(A) are listed on table 2, and 
/i(A) = 1 — 6(A) + c?i(A). In ref. [|] such vertex model was mapped on a certain isotropic 
branch of the Izergin-Korepin model, and consequently its algebraic Bethe Ansatz solution is 
similar to that developed by Tarasov |27| ( see also ref. |J ). By including the inhomogeneities 
{Ai, • • • , A n } we conclude that the corresponding eigenvalues are 
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where the variables {/ij} satisfy the equation 
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Now we come to the point where all the results can be gathered together in the following 
way. Supposing we are interested in the eigenvalues of our original vertex model, we start 
with the eigenvalue formula (55) and use recurrence relation (65) until we reach the problem 
of diagonalizing the Bi, C±, D 2 and Osp(l\2) models. Then, we have to take into account our 
results for the eigenvalues of these systems, which are collected in equations (68,71,75,78). By 
using this recipe, it is straightforward to find the eigenvalues expressions for the vertex models 
we have so far discussed in this paper. We now will list our results for the eigenvalues. For the 
B„ (n > 1) model we have 
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For C n (n > 2) we have 
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where the functions G t (X, {A^}) are 
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For D n {n > 3) we have 
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where the functions G t (X, {A^ }) are 
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where the functions G/(A, {A^}) are 
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Finally, for Osp(l\2n) (n > 1) we have 
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We see that the eigenvalues depend on the parameters |{A^}, {A^}, • • • , {A^ 1 }}, which 
represent the multi-particle Hilbert space of the many steps necessary for the diagonalization 
of the nesting problem. In the expressions (80-93) we have already performed convenient shifts 
in these parameters, {A^} — > {A^} — in order to present the final results in a more 
symmetrical way. In table 4 we show the values for the shifts These set of variables are 
constrained by the Bethe Ansatz equation, again at each level of the nesting. The same recipe 
described above for the eigenvalues also works for determining the corresponding Bethe Ansatz 
equations. We start with equation (57) and each step of the nesting is disentangled by using 
the recurrence relation (66). When we reach the last step, we take into account the Bethe 
Ansatz results (69,72,76,79) for the inhomogeneous Bi, Ci, D 2 and Osp(l\2) vertex models. It 
turns out that the nested Bethe Ansatz equations have the following structure. All the vertex 
models share a common part, which can be written as 
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where s(n) = n — 1 for the B n , Osp(2n — 1|2) and Osp(l\2n) models ; s(n) = n — 2 for the 
C n , Osp(2\2n — 2) models; s(n) = n — 3 for the D n and Osp(2n — 2|2) models. The remaining 
equations are somewhat model dependent and below we list their particular forms. For the 
models B n (n > 2), C n (n > 3), D n (n > 4) and Osp(l\2n) (n > 2) the equation for the first 
root {A^} is given by 



\f - 1/2 
X? + 1/2 



n 



A 



(i) 



xf -W + ik" xy> - xr - 1/2 



1 m 2 x(2) 
1 -Q A j ~ A k 



+ 1/2 



(i) 



(2) 



(95) 



while for the Osp(2n - 1|2) (n > 2) , Osp(2\2n - 2) (n > 3) and Osp(2n - 2|2) (n > 4) we 
have 
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Due to the peculiar root structure of the models C 2 , D 3 and Osp(2\2), Osp(A\2), their first 
Bethe Ansatz equations are a bit different than that presented in equations (95) and (96), 
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respectively. In order to avoid further confusion we have to quote them separately. For the C2 
model the Bethe Ansatz equation for {A^ } is 



A 



(i) 



1/2 



\ { - ] + 1/2 



mi 

n 



X ( P - At Xj - 1 



(i) 



k=i, k^j Aj 



(i) 



rn 2 

n 



A 



(i) 



Af + 1 



J_ 1 11 \W 



A 



(2) 



(97) 



while for the Osp(2\2) model we have 
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Furthermore, the first equation for the D 3 model is 
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and for the model Osp(A\2) we have 
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The Bethe Ansatz equation for the last variables are common for the pairs {B n , Osp(2n - 
1|2)}, {C n ,Osp{2\2n-2)}, {D n , Osp{2n-2\2)} and are given as follows. For B n and Osp{2n- 
1|2) models the variables {A^ } satisfy the equation 
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For C n and Osp(2\2n — 2), the parameters {A^ n ^} and {A^ } satisfy the equations 
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For n = 2 ( C2 and Osp(2\2) models ) we have to consider only the last equation given in (102) 



30 



For D n and Osp(2n— 2|2), the parameters {A^ n 2 ^}, {A^ } and {A^- •*} satisfy the equations^] 
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and finally for the Osp(l\2n) vertex model the equation for variables {A^ } isQ 
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We would like to close this section with the following remarks. We begin by discussing 
and comparing our results to previous work in the literature. It is not difficult to check that 
our results for the Bethe Ansatz equations are equivalent to the analyticity of the eigenvalues 
as a function of variables {A^}, • ■ ■ , {A^- }|, i.e. all the residues on the direct and crossed 
poles vanish. This is precisely the main important ingredient entering in the framework of the 
analytical Bethe Ansatz. Indeed, the results of this paper for the eigenvalues and the Bethe 
equations of the B n , C n , D n and Osp(l\2n) models are in accordance to those obtained from 
the analytical Bethe Ansatz approach in refs.0, |3|, f|, ||. As we have already commented in the 



introduction, the D n vertex model has also been solved by the algebraic Bethe Ansatz in ref.||14j|. 
These authors argued that the nesting problem, in a convenient basis, could be transformed to 
two commuting eigenvalue problems. One of them has the permutation operator as the main 
intertwiner, and the algebraic Bethe solution goes along to that known to work for the multi- 



state 6- vertex models [[11], [12], y| . The other one was related to the Temperley-Lieb operator ( 
in our notation of section 2 ), no explicit algebraic solution was attempted, and the eigenvalue 

10 Here we have assumed the identifications A^™ ^ = A^ and A^™' = A^ Moreover, for n = 3 (D3 and 

Osp(4\2) models ) only the last equation in (103) matters. 

n We remark that the Bethe Ansatz equations and eigenvalues for the models B\, C±, D2 and Osp(l\2) can 

be obtained from equations (68-69,71-72,75-76,78-79) by setting the inhomogeneities to zero. 

12 We remark that ref.j|] have used the grading /•••/&/•••/ and here we have taken the sequence 

b ■ ■ -bfb ■ ■ -b. This is the reason why the phase factors of ref.jfj] and that of equation (104) are different. 

We also have noticed misprints in rcf. |2] concerning some eigenvalue expressions. 
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results were obtained via the crossing property. In our approach, however, we deal with these 
two operators together in the diagonalization problem and the explicit use of crossing it is not 
needed. We also noticed that while our formulation works for vertex models having both even 



or odd numbers of states q per link, the basis used in ref. [14| appears to be suitable only for q 



even. Since our approach and the one of ref. jnj are quite different from the very beginning 



we were not able to find a simple way of comparing the results for the eigenvectors. 

For the models based on the superalgebras Osp(n\2m), we notice that the results for the 
eigenvalues and Bethe Ansatz equations present some additional phase factors when compared 
to that obtained for the pure "bosonic" B n , C n and D n models. These phases are the only liquid 
difference between the standard and the graded formulations of the quantum inverse scattering 
method. They are indeed twisted boundary conditions (periodic and antiperiodic) having a 
fermionic sector dependence ( see e.g. ref. Q )• In appendix C we show how such phase 
factors can be absorbed in the supersymmetric formulation of the Osp(2n — l\2),Osp(2\2n — 
2), Osp{2n — 2 1 2) and Osp{l\2n) vertex models. Without these extra factors, the nested 
Bethe Ansatz equations becomes more symmetrical, similar to what happens for the "bosonic" 
models. It turns out that such symmetrical Bethe Ansatz equations can even be formulated 
in a more compact form, reflecting the underlying group symmetry of these vertex models. In 
order to see this, we first need to make the following definitions. Let C a b be Cartan matrix 
associated to the Dynkin diagrams of figure 1, and rj a the normalized length of the root of type 
/3 a . Here we assume that the length of the long root is 2. This means that r\ a = j- 2 a > and 
rj a — 1 for a long root []]. Taking into account these definitions, it is possible to rewrite the 
Bethe Ansatz equations (94-103) as0 



\(a) _ <5q,i ' 

j 2r,a 

\[a) , $a,l 

X 3 +2^T. 



r(n) m b ^( a ) _ ^( b ) _ C a,b 

II II X) ^W~Sj > 3 = l,...,ma) a = l,...,r(n) (105) 



6=1 k =i, k^j Aj - A k + 



13 More precisely, for D n and Osp(2n — 2|2) r\ a — 1, a = l,...,n; for B n and Osp(2n — 1|2) r\ a = 1, 
a = 1, . . . , n — 1, and rj n = 2; for C n and Osp(2\2n — 2) r\ a = 2, a = 1, . . . , n — 1, and r\ n — 1 (see for instance 
ref.f§). 

14 We notice that for the C n and Osp(2\2n — 2) vertex models one should rescale all {A^} by a factor 2, 

(a) {A^ a ^} 

{Xj } — > — 7j — in equation (105) to recover the previous results (94,95,97,98,100). 
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where in general r(n) is the number of roots of the underlying algebra. For the vertex models 
solved in this papers r{n) = n. 

The idea that the Bethe Ansatz equations can be compactly written in terms of their 
corresponding Lie algebra goes back to the work of Ogievestsky, Reshetikhin and Wiegmann 
[|30|j who have conjectured similar formulas for the Bethe Ansatz equations of factorizable S- 
matrices based on the standard Lie algebras A n , B n , C n , D n , Eq, E 7 , E$, F 4 , G2 ■ Thus, 
our results can be seen as a rigorous proof of this conjecture for the rational B n , C n , D n 
vertex models. Moreover, together with the results of ref. || for the Sl(n\m) algebra (see also 
ref. |29|), they also show that this conjecture can be extended to the case of superalgebras. An 



exception in this construction is again the Osp(l\2n) vertex model. In order to fit equation 
(105), we need to give a special meaning to the black bullet • of figure 1 ( root {A^} in the 
Bethe Ansatz ). One could interpret it as a peculiar two-body self interaction for root {X^}, 
as the one present in the right hand side of equation (104). 

Our last remark is concerned with the underlying spin chain associated to the braid-monoid 
vertex models. The i?-matrix (4) is regular at A = 0, thus the local conserved charges can 
be obtained through the logarithmic expansion of the corresponding transfer-matrix around 
A = 0. The first charge is the momentum itself (we are assuming periodic boundary conditions) 
and the next one, i.e. the logarithmic derivative of T(A) at A = 0, is the Hamiltonian. From 
the expression of the i?-matrix -R(A) we found 

'H = j2P? + ]rE l (106) 
i=i lA 

Analogously, the eigenenergies E(L) of the corresponding Hamiltonian (106) can be calcu- 
lated by taking the logarithmic derivative of A(A, {A^}, . . . , {A,- }) at the regular point A = 0. 



The eigenenergies E(L) are parametrized in terms of the variables {\\ } byQ] 

mi 1 

E(L) = ±V^ ±L (107) 

h [A, (1) ] 2 - 1/4 



15 It has become a tradition in the literature to normalize the Bethe Ansatz variables by a pure imaginary 

Til -, 

r - 

^[A 4 (1) ] 2 + l/4 



\W 7711 ] 
factor as — > -4—. In this case, equation (107) reads E(L) = T J]] — ^pj ± L 
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where the signs ± is related to the two possibilities for the Boltzmann weights of type a(A), 
i.e., a(A) = 1 ± A. The variables {\\ } satisfy the general nested Bethe Ansatz equations 
(94-104), and therefore the spectrum of all Osp(n\2m) chains can be determined solving these 
equations. A preliminary study of the root structures of the Bethe Ansatz equations (94-104) 
have shown an intricate behaviour for the general Osp(n\2m). We leave a detailed analysis of 
these root structure, as well the ground state and the low-lying excitations for a forthcoming 
paper. 

6 Conclusions 

In this paper we have solved exactly a series of rational vertex models based on the braid- 
monoid algebra from a rather unified perspective. The general construction for the i?-matrices 
presented in section 2 played an important role in many steps of our formulation of the quantum 
inverse scattering method. In particular, it is notable how such construction becomes useful to 
determine the universal character of the fundamental commutations rules presented in section 
3. It appears natural to us to blame the properties of the braid-monoid algebra as the main 
mathematical structure behind such general picture. 

The commutation relations allowed us to determine the eigenvectors and the eigenvalues 
of the transfer matrix of the vertex models invariant by the B n , C n , D n , Osp(2n — 1|2), 
Osp(2\2n — 2), Osp(2n — 2|2) and Osp(l\2n) symmetries from systematic point of view. As a 
consequence, we have been able to establish the solution of the general Osp(n\2m) invariant 
spin chain. Their Bethe Ansatz equations have been formulated in terms of basic properties 
of the underlying group symmetry such as the root structure. This means that, in principle, 
by solving equations (105) together with (107) we can determine the spectrum of the general 
Osp(n\2m) spin chain (106). 

One possible extension of this work is to consider the trigonometric analogs of the vertex 
models solved in this paper. In general, we expect that the trigonometric vertex models based 
on the Birman-Wenzl-Murakami algebra [|IJ |2(| can be solved by introducing few adaptations 
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to our construction of sections 3 and 4. In addition, we have noticed that similar ideas also 
works for the one- dimensional Hubbard model . In this case, however, some other peculiarities 
such as the non-additive property of the i?-matrix need to be disentangled, too. Due to the 
special character of the Hubbard i?-matrix and its recognized importance in condensed matter 
physics, we have dedicated a brief account of our results in a separated publication [pB| . It 
remains to be seen whether these are isolated examples or our construction could even be more 
widely applicable that one would think at first sight. The last possibility would suggest that 
a more deep and essential mathematical structure is still to be grasped from the Yang-Baxter 
algebra. 
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Appendix A : The two-particle state 

In this appendix we present some details concerning the complete solution of the eigenvalue 
problem for the two-particle state. Here we shall need few extra commutation relations between 
the fields 5(A), B*(X), C(X) and C*(X). They are given by 

C a (\)B b (n) = B b (fi)C a {\) - T7T- — c [B(X)A ab (fi) - B{^)A ah {X)} (A.l) 
B* a (X)B b (») = B b {n)B:{\) - —±—[F{\)A ab {ii) - F(/i)A o6 (A)] (A.2) 

C:{X)BM |0> = U ^Zrf PO^A) ~ A aaWAM] |0) (A.3) 

where a,b = 1, • • • , q — 2. 

We start our discussion by considering the wanted terms for the two-particle state. Since 
the terms proportional to Bi(\^)Bj(\^)J r: ' 1 have been already cataloged in section 4 ( see 
expression (55) ) we turn our attention to those proportional to — ^^ F^X^^B^X^)^.^ 
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. In other words, we have to show that the wanted terms proportional to the second part of 
the two-particle eigenvector (40) produce the eigenvalue expression (55) for mi = 2. It is 
not difficult to collect the contributions proportional to [a(A)] L and [e„(A) L ]. They are given 
respectively by 

" a(Ai 1} - A) 1 a (Aj 1 )-A)rf n (AS 1) -A)e n (AS 1) -A( 1) ) 

_e n (A< 1) - A) 6(A?> - A) b(X? - A) e^A^ - A) d n (X? - A?>) 

" a(A - X?) _ 1 d n (X-X^)e n (X^ -X ( i ] ) 
e n (X - X?) e n (X - X?) e n (X - X?) d n (X? - A? } ) 
Now if we use that the Boltzmann weights satisfy the following identities 



[a(X)] L (A.4) 



[e n {X)\ L (A.5) 



(i) 



a(X? - A) + 1 a (\V-\)d n (XY>-X)e n (XY>-X' 2 1> ) 
e n (X? - A) 6(A? } - A) b(X? - A) e^A^ - A) d n (X? - A^) 



1=1 



a(X-X\ 



(l)^ 



1 ^(A-A^K^-AH 



_fr KA-A t (1) ) 
tie^A-Af)) 



.e„(A-A«) e n (X-XP)e n (X-XP)d n (xU-XP)_ 11 " ^ 
it is immediately seen that these two contributions are precisely the first and the second terms 
of equation (55) for m\ = 2. 

The analysis of the wanted term proportional to [&(A)] L is a bit more involved. The basic rea- 
son for that is the presence of some "asymmetric" terms of kind [b(\)a(\2^)] L F(\i^) Y^ij^ijF^ 
and the peculiar form of the nested eigenvalue A^(A, {X^}). We find that the terms contribut- 
ing to this last wanted piece are given by 



"(9-2) 



1 - 



[b(X - X[ 



(1)M2 



(1) 



X 



[b(X)a(X^Y^ 



A?) 



e n(Ai — A 2 



V) 



[b(X)a(X^)\ 



W\1L 



d n (X-X?) 



[b(X - a[ 1} )] 2 e n (A - A^) tjhn 



Eu-[^ (1) (A-AS 1) )]irF(Ai 1) )^, 



(A.8) 

(A.9) 
(A.10) 



e n (A-A«) b(X-X^) b(X-X^) 
where the first term (A.8) comes from the eigenvector part proportional to F(X^) while the 
remaining ones originate from the term B i (X^)Bj(X^)J : ' jl - In order to go on we need to take 
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advantage of the identity 



A«(A, {X?})^ 1 = [T«(A, {A^})]^- 



(A.ll) 



By expanding expressions (A.8-A.10) with the helping of identity (A.ll) we have checked 
that these three terms together produce the third piece of the eigenvalue expression (55). 

We now turn our attention to the unwanted terms. They are generated when the parameter 
A of the diagonal operators B(X), A aa (X) and D(X) exchanges with the variables {A^A^} 
parametrizing the eigenvector $ 2 (Ai 1 ' ) , A^)^. Basically, we have three kinds of such terms and 
they are proportional to 

B^B^), BXXIB^), F(X) (A.12) 

When A- 1 - 1 = A^ we just have two contributions for the first two unwanted terms (A.12), 
and they are canceled out by same procedure presented in section 4 ( see equations (51,52) ). 
By contrast, when A- 1 "* = A^, several terms appear and the simplifications are more involved. 
For example, in the case of the unwanted term B i (X)Bj(X^)J 7: ' 1 ' we have four contributions 
coming from S i (A[ 1) )S i (A? ) ) which are given by 



( j hl a(A - A^ 1} ) 



6(A - X?) b(X - \ { 2 L> ) % 



(i) 



[a(A^)] 



i 



b(X - aS 1} ) 6(aS 1} - A^ 1; ) 



6(A? ) - A) b(X^ - A) t^ra 
+[a(A?>)f - — 



1 



^B^B^)^ (A.13) 



6(A« - A) 6(A?> - A?>) % 
and only one coming from the term F(A 1 1 ' ) )^*.JF, which is given by 



[a(A?>)] 



l)ML^n(Ai^ — A) C? n (Ai ±J — A^) 



(1) 



e n(X[ ^ ~ A) e n (X\ ' — X^') ijlm 

The simplest way of checking that the terms (A.13) and (A. 14) together cancel out is by using 
explicitly the Bethe ansatz equations (45) for variable A^, namely 

1 



(A.14) 



(i) 



(i) 



i(4 1) )] L E[^ (1) (Af ) -A? ) )]^ i (A.15) 
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Substituting this expression in the first two terms of equation (A. 13), all the unwanted 
terms become proportional to [a(A2^)] L . Then, we have verified that such final expression gives 
a null result with the helping of some additional identities between the Boltzmann weigths. 
Analagously, a similar procedure can be used to deal with the unwanted term B*(X)Bj(X^). 
The five contributions for this unwanted term are collected below 

1 d n (X - X?) 



[a(\?>)] 



b(X — X[ ^) e n (A — A 2 ') aijlm 



E t^ (1) (A - Ai 1) )]2?s:(A)S I (Ai 1, )^ i ; 



e n (X - XI 1 ') &(A?> - Ai 



-[&(A? } )] L 
-[a(A«)] L 



6(A-Aj 1) ) 4(A-A^) rf n (A-Ai iJ ) 



e n (A - A?>) e n (A - A?>) e n (X - A«) 6(A[ 1J - A«) 
^E^AWAiV^ 



1 cj ra (Aj 1} - A ?>) 
6(A-A( 1) )e„(Ai 1) -A 2 j ijt 



^B* l {X)B l {X i l ) )^\ 

ijl 

(A.16) 



Lastly, we have nine contributions to the third term proportional to F(X). The first seven 



come from B i (\i'^)Bj(X2' > ) and the other two from F(X\ L> ). An easy way of verifying that all 
these nine terms cancel out is by using explicitly the Bethe ansatz equations (45) both for A^ 
and X^ ■ The final expression can be disentangled only in terms of [a(X ( [ S) )a(X^)] L F(X), and 
we have checked that it gives a null by using Mathematica™ . 

Appendix B : The three-particle state symmetrization 
This appendix is concerned with the symmetric properties of the three-particle state. We 



0-h 



begin our discussion with the A^ <-> X^ 1 permutation. After this permutation the (vector) 
three-particle state (46) looks like 

l3(A?\Ai 1 \4 1 >) = S(A«)®$ 2 (Af),^) + 

[^F(X?)B(X?)B(X?)} MAMU 1 )) + [f® F(A«)5(A«)5(A«)] h 2 (X?\ A«, A?>) 

(B.l) 



This vector can be related to ^(A^, A 2 *\ A^) by commuting the fields B^ 1 ') and B(X 
in expression (46). In order to do that, we have to use the commutation rule (29), which now 



0-h 
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reads 



a(A 2 — A3 ) 



^n(A 2 ^ — A3 "*) 



F(A«)B(A«) - {1 ~^ W ~f )} F(A.W)B(A?») 
a(A 2 - A 3 ; ) 



(B.2) 



By substituting this relation in equation (46) and by comparing the final result with the 
vector (B.l), we find that the symmetric property (50) holds for variables {X^\ A 3 ^} provided 
the functions h^X? , X? , \?) and h 2 {X^ , X^ , X^) satisfy 



h 2 (X?\ A?\ A?>) = MAj 1 ), A^A^ p'Yn /n't' (B.3) 



A 23 l A 2 ~ A 3 J 

a(A?> - A?>) 



MM > a 2 , A3) = ^(a^, ^W rXm SI (B-4) 



,(1) \(i)a. _ L/iW \W a(i)a A 2 3 (^2 -A 3 

a(A^ 1} - A3 

In fact, the two expressions above are equivalent due to the unitarity property of the 
auxiliary matrix X^(X), namely 

-u) = a(u)a(-u)I (B.5) 

Therefore, in order to determine these functions we still need to find an extra constraint. 
This comes out if we perform the permutation between variables A^ and A 2 *\ The procedure 
is similar to the one described above, but now some additional commutation rules are needed. 
For example, besides commuting the creation fields B(X^) and ^(A^), we have to commute 
S(aS 1} ) and F(A^ 1} ) and also F(a! 1} ) and B(X^). This latter step is worked out with the 
helping of commutations rules (31) and (32), respectively. We then are able to determine the 
function ^(A^A^A^) by imposing the necessary restriction that the terms proportional 
to |*<g) F(X^)B(X^)B(X I 2^) need to be canceled out. We find that this condition is verified 
provided 

i aW _ "nlA - X 2 ) a(A 3 - X 2 ) 

tii^ , A 2 , A 3 ) — ttt rm7T7T7ri TiTu^ ^-"J 



and by substituting (B.6) in (B.3) we are able to fix h 2 (Xi'\ X^\ A^ 1 - 1 ) ( see equation (48 
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Besides that, some other non-trivial checks need to be performed in order to verify the 
consistency of the procedure mentioned above. In the process of using the three commutation 
rules mentioned above, certain extra identities emerge, and they are given by 



[£®B()g ) )F(\Y>)B(\)t>)]b(\)t>->b 



>(i) 



e n lA 2 — A 3 ) 



(1) 






2 


) 






d n 


(A 


(i) 
i 


-A«) 




(A 


(i) 
l 





jn(A^~4_) , .,!!..) ,(11 .UK 

+ "2 (A >^2 ; A 3 J 



-A^) e„(A^-A^0 



(i) x (i) x (ih 

1 ) A 2 ) A 3 t 

(B.7) 



and 



-6(A« 



cj„(A^ — A^ ; ) 



KA«-A«) e„(A«-A«) <»-AW) e^-A™) 



^ ) )f s ^^[^V(Af ) )S(Af ) )®a = 



\i ® F{\f)B{\V)B{\V)\ [-h 2 (\\\^>, + - WfoW, \ ( i>, A 3 ±J )J (B.8) 

In order to show that these relations are satisfied, we end up proving remarkable properties 
between £ and the auxiliary matrix X^\u). They are listed below as 

dJu) e n (-u) 



,(1) 



x (i) 



[6 2 5 3 (y)]x«( M ) 
[5i(y)6 3 ]xgV) 



d n (u) e n (-u) 



e n (u) d n (-u) 



[l-tb(u)][S 12 B 3 

[i-Hu)]^)^ 



u 



u - AW 
u 



[£iA{y)\ 



[Z 12 B 3 (y)]X£(u) = [^ 12 B 3 (y)]+u[^ 2 B 3 (y)P^ 9 ] ~ ^3^iy [Bi(y)^\ 
[£ 12 B 3 (y)Pg9]X$(u) = [i 12 My)P& )9 ]+^Bi(y)^3] - ^z^i)*[fi2Ss(y)] 



(B.9) 

(B.10) 
(B.ll) 
(B.12) 
(B.13) 



where the lower indices indicate the position where vector £*, vector B(y) and matrix X^(u) 
acts in a non-trivial way. P^ 9 means the permutator entering in the construction of JW(A). 
Besides that, we also note the following useful identity £i 2 B 3 (y)P 23 g = Bi(y)£ 23 P]^ 9 '. 
Appendix C : The graded quantum inverse approach 
The purpose of this appendix is to discuss the main modifications occurring in the commuta- 
tion rules, eigenvalues and Bethe ansatz equations for the models Osp(2n — 1|2), Osp(2\2n — 2), 
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Osp(2n — 2 1 2) and Osp(l\2n) when we formulate their solution in terms of the graded inverse 
scattering framework |TTj|. In this formalism, the vertex operator Cj^i(X) entering in the mon- 
odromy matrix (12) satisfies the graded Yang-Baxter equation and is given by acting the graded 
permutation operator P 9 on the i?-matrix (4), namely 

Z#(A) = (-l)^ 6 >i^(A) (C.l) 

In order to accomplish this change, the integrability condition (14) now reads 

R{\, ii)T(X) <§> T(/i) = T(/i) <§> T(X)R(X, /i). (C.2) 

s / s \ cd 

where the symbol ® stands for the supertensor product I A ® B ) = (— l) p ^^ p ^ +p ^ AacB^. 

\ J ab 

Moreover, the transfer matrix T(A) is written in terms of a supertrace of the monodromy 
matrix T(A) as 

T(A) = StrT(X) = 2(-l) p(a) ^a(A) (C.3) 

a 

The modifications (C.l— C.3) are responsible by the appearance of extra signs in some terms 
of the commutation rules, in the "diagonal" conditions (19) and in the eigenvalue problem (18). 
Furthermore, the way that such signs enter on these relations depend much on the grading 
sequence we have chosen. We recall that we have taken the fb- ■ - bf grading for the first three 
models and b- ■ -bfb- ■ -b for the Osp(l\2n) vertex model. Therefore, we have to present the 
modifications for these two classes of grading separately. We begin by listing the results for the 
Osp(2n — X 1 2) , Osp(2\2n — 2) and Osp(2n — 2 1 2) models. Concerning commutation rules, the 
extra signs appear only in relations (22,23,24,29) and for the following right hand side pairs of 
terms 

BiBj -> -BiBj, B, t B -> -BiB, B, t D -> -B, t D, B*B -> -B*B. (C.4) 

All the other fundamental commutation relations remain unchanged. Furthermore, the 
action of the diagonal operators B(X), A ab (X) and D(X) on the reference state (15) becomes 

B(X) |0> = [-a(X)] L |0> , D(X) |0> = [-e n (X)] L |0> , A aa (A) |0> = [b(X)] L |0> , a = 1, . . . , q - 2. 

(C.5) 
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while the eigenvalue problem (18) becomes 

l-B(\) + Y J A aa (\) - D{\)\ |$) = A(A) |$) (C.6) 
a=i 

Taking into account these relations one can easily verify that the phase factor (— i) L -" t i- 1 
present in the Bethe ansatz equations of these models is canceled out. For sake of completeness 
we also present the final results for the eigenvalues in the graded formalism. For Osp(2n — 1|2) 
(n > 2) we have 



A Os P (2n-l\2) (A;{Af},...,{A5 n) } 



mi \ _ \ (!) I 1 mi \ _ \ (!) _i_ _ q 2n-l 

j=i A — A — 2 j=i A — A • + n — z ; = i 



(C.7) 



and for Osp(2|2n — 2), (n > 2) we find 



11 A - A (1) + - 

=1 A — Aj - 2 j=l A — Aj + -j- ,=i 



mi \ _ \(!) , 1 mi \ _ _|_ 2n-3 2n-2 



(C. 

while for Osp(2n — 2|2) (n > 3) we have 



A S p(2n-2|2) ( A;{A (. 1 ) }) ... ){A W } 



(!) , 1 mi \ \W i 2n-7 2n-2 



mi \ _i_ l mi \ _ _i 

-i-°w\ L n , % i - i-^\)] L n %; 2 ^ 5 + ika)] l e <?,(a, {a?> » 

j=l A — Aj — 2 j=i A — Aj H 2^ — 1=1 

(C.9) 

where functions G/(A, {A^ 1 ) remain unchanged and are given by equations (87), (89) and (91), 
respectively. 

The analysis for Osp(l\2n) model is a bit simpler because we just have one "central" 
fermionic species in the grading sequence. In this case, we have only extra signs for the terms 
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BiAji and [B*]iAji in relations (22) and (24), respectively. The change in the eigenvalue (92) 
now appears on the n-th term, and the final result is 



A o sp( i|2„) ( A; {A< 1 >},...,{A<">}) 



i"W) L ff \ %J + MA),i n Vvt" +1 + 

j=i A - A] + ± j= i A - A) + n 

n-1 2n-l 

[b(\)] L £ Gi(A, {Af }) + [b(X)] L £ G,(A, {Af }) - [6(A)] L a t (A, {Af } }) (CIO) 

2=1 2=n+l 

where functions G/(A, {Af } for Z = 1, • • • , n — 1, n + 1, • • • , 2n — 1 are still given by (93), but 
G n (A,{Af}) is 

m (n) ( , _ (n-lK ^ _ (n+2) v 

G»(A, {Af }) = n ' J„ + .' )(A ^ 2 ' (Cll) 

Analogously, the Bethe ansatz equations only modify for the last root {Af }. Instead of 
equation (104) we now have 

"fx 1 - Aj"" 1} + 1/2 (Af - Af + l/2)(Af - Af - 1) _ i 
A=i - A^ - 1/2 M (Af - Af - l/2)(Af - Af + 1) 
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Tables 



Table 1 : The number of possible states per link q and the parameters t and JC for the 
models B n , C n , D n , Osp(2n - 1|2), Osp{2\2n - 2), Osp{2n - 2|2) and Osp(l\2n). 







C n 




Osp(2n- 1|2) 


Osp(2|2n-2) 


Osp(2n-2|2) 


Osp(l\2n) 


Q 


2n+l 


2n 


2n 


2n+ 1 


2n 


2n 


2n+ 1 


K 


2n+l 


-2n 


2n 


2n-3 


4-2n 


2n-4 


l-2n 


i 


1 


-1 


1 


1 


-1 


1 


-1 



Table 2 : The main five Boltzmann weights for the vertex models B n , C n , D n , Osp(2n — 1|2), 
Osp(2\2n - 2),Osp(2n - 2|2) and Osp(l\2n). 







Cn 




Osp(2n- 1|2) 


Osp{2\2n-2) 


Osp(2n -2|2) 


Osp{l\2n) 


a{\) 


A + l 


A + l 


A + l 


1-A 


1-A 


1-A 


A + l 


6(A) 


A 


A 


A 


A 


A 


A 


A 


c„(A) 


n-l/2 


2A+n+l 


n-l 


2A+n-5/2 


n-l 


2A+n-3 


2A+n+l/2 


A+n-1/2 


A+n+1 


A+n-1 


A+n-5/2 


A+n-1 


A+n-3 


A+n+1/2 


4(A) 


—A 


A 


A 


-A 


A 


-A 


A 


A+n-1/2 


A+n+1 


A+n-1 


A+n-5/2 


A+n-1 


A+n-3 


A+n+1/2 


e„(A) 


A(A+n-3/2) 
A+n-1/2 


A(A+n) 
A+n+1 


A(A+n-2) 
A+n-1 


-A(A+n-3/2) 
A+n-5/2 


A(A+n) 
A+n-1 


-A(A+n-2) 
A+n-3 


A(A+n-l/2) 
A+n+1/2 
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Table 3 : The main five Boltzmann weights of 'nested' matrix X^ k \ They are the same for 
the pairs {B n , Osp(2n-l\2)}, {C n , Osp(2\2n-2)} and {D n , Osp{2n-2\2)}. The corresponding 
crossing parameter A^ fc ^ is also listed. 



I«(A) 


{B n ,Osp(2n -1\2)} 


{C n ,Osp(2|2n-2)} 


{D n ,Osp(2n -2|2)} 


Osp(l|2n) 


a(A) 


A + l 


A + l 


A + l 


A + l 


6(A) 


A 


A 


A 


A 


Cn-fe(A) 


n-fc-1/2 


2A+n-fc+l 


n— fc— 1 


2A+n-fc+l/2 


A+ra-fc-1/2 


A+n-fc+1 


A+n-fc-1 


A+n-fc+1/2 


^n-fc(A) 


—A 


A 


-A 


A 


A+n-fc-1/2 


A+n-fc+1 


A+n-fc-1 


A+n-fc+1/2 


e„-fc(A) 


A(A+n-fc-3/2) 
A+n-fc-1/2 


A(A+n-fc) 
A+n-fc+1 


A(A+n-fc-2) 
A+n-fc-1 


A(A+n-fc-l/2) 
A+n-fc+1/2 




-ra + k+ 1/2 


— n + fc — 1 


-n + fc + 1 


-n + k- 1/2 



Table 4 : The shifts 5^ performed on the Bethe Ansatz variables A^ . 





B n and Osp{l\2n) 








(3/2, (3 = l,---,n 


1 


■(5/2, (3 = l,---,n-l 
. (n + l)/2, /3 = n 


! 


'/3/2, /3 = l,.-.,n-2 
. (n-l)/2, /3 = ± 





Osp(2n- 1|2) 


Osp(2|2n-2) 


Osp(2n-2|2) 




(/3-2)/2, /3 = 1,.-., n 


! 


'(/3-2)/2, /3 = l,.-.,n-l 
. (n-l)/2, = n 


! 


'(/3-2)/2, /3 = l,-..,n-2 
. (ra-3)/2, /3 = ± 
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Figure 

Fig. 1. The Dynkin diagrams of the Lie algebras B n , C n and D n and of the superalgebra 
0sp(n\2m) (e.g see ref. p5) ). The symbols 0> ® an d • stand for the simple roots of the 
algebras 5/(2), 5/(111) and Osp(l\2), respectively. 
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